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Abstract : An SU(3) symmetric model with high predictivity for octet meson 
(tt, K) quark fragmentation functions with a simple approach to SU(3) symmetry 
breaking (due to the relatively heavy strange quarks) is extended to the singlet sector, 
with some reasonable assumptions, in order to study 77 and 7/ fragmentation. Due 
to the presence of SU(3) symmetry, fits to the tt and K data help to constrain the 
fragmentation functions for 7/ and rf mesons. The use of 2-jet and 3-jet (especially the 
gluon jet) inclusive meson production in e + e~ collisions and tt, rj inclusive production 
in pp collisions enable the extraction of the gluon fragmentation functions as well. 
While sea quarks in strange mesons (K ± , K°, K ) and the heavier 77, rj mesons 
are suppressed by a factor of \h = m^/m 2 H < 0.1 for H = K,r],r)', the gluons are 
not as severely suppressed: fjf ~ 0.3-0.35. A detailed parametrisation of the three 
independent fragmentation functions V(x,Q 2 ), j(x, Q 2 ), and D g (x,Q 2 ) are given in 
>■ ■ LO QCD; all nonet meson gluon fragmentation functions and quark fragmentation 
functions of all flavours can be expressed in terms of these functions with a few model 
! parameters including A# and f g . The data prefer a nonet mixing angle —24° < 6p < 
— 16° and rule out no-mixing. 

3 ' PACS Nos: 13.87.Fh, 13.65. +i, 13.60.Le, 14.40.Aq, 13.85.Ni 

o\ : 
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> ! 1 Introduction 

X ' 

The PHENIX experiment at RHIC [lj has measured inclusive transverse momentum 
spectra of rj mesons in the range pt ~ 2-11 GeV at mid-rapidity (\y\ ~ \rj\ < 0.35) 
in pp, d Au and AuAu collisions at ^/snn = 200 GeV. Combining with earlier data 
on tt production [2j, as well as other similar data [3], their analysis shows that the 
77/71-° production ratio is roughly constant, ~ 0.5, for a range of center-of-mass energies 
(y/$NN ~ 3-1800 GeV). This number agrees with e + e~ annihilation data at y/s = 91.2 
GeV for high scaled momentum x p . This ratio is interesting because of its potential 
as a good signal for quark-gluon-plasma (QGP). Hence, a good base-line study of r\ 
fragmentation functions in pp collisions is essential for use as a normalisation factor 
in the study of this signal of QGP. 

There are many previous studies on meson fragmentation functions at both leading 
order (LO) and next-to-leading order (NLO) in QCD. While the fits at NLO on tt and 
K inclusive production presented in Ref. |4j are classic and detailed, several updated 
approaches have been studied [5]. While these have concentrated on the dominantly 
produced tt and K mesons which constitute the bulk of the particle content of jets, 
for example, about 91% at the Z-pole, 7/ (and 7/) fragmentation has not been much 
studied, as has been pointed out in Ref. [I]: one such study is Ref. [6]. 
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In our earlier work [7], an SU(3) model of light quark (u,d,s) fragmentation 
functions for octet baryons and pseudo-scalar mesons which is in good agreement with 
hadro production (tc and K) data in e + e~ annihilation experiments was developed. 
Due to the presence of SU(3) symmetry, the model can predict the (purely octet part 
of the) r) fragmentation functions; in fact, the pure octet rjg fragmentation functions 
were listed in this paper. However, due to octet-singlet mixing, the r\ and rj' mesons 
were not considered in this paper. The RHIC/PHENIX result motivates us to extend 
this model to the singlet sector of the pseudo-scalar mesons. 

One of the main advantages of our approach is the presence of SU(3) symmetry 
so that the valence and sea quark fragmentation functions of the octet mesons are 
related to each other and can be written in terms of just three independent fragmen- 
tation functions, a valence function V(x, Q 2 ), a sea function j(x, Q 2 ) and the gluon 
D g (x,Q 2 ). The model thus has high predictive power. SU(3) breaking effects and 
singlet-octet mixing are then included by introducing several parameters, which are 
taken to be constants in this simple-minded approach and yet give good agreement 
with data. In particular, it is possible to greatly improve the quality of the fits by 
using the abundant data on 7r and K inclusive production to fix the quark and gluon 
fragmentation functions and so to obtain both the gluon fragmentation function and 
quark fragmentation functions for all flavours for the f] meson. 

In the next section, we briefly review the cross-sections and rates for inclusive 
particle production in e + e~ (2-jet and 3-jet) and pp collisions. In section 3, we 
highlight the salient features of the model, which is based on SU(3) flavour symmetry 
of the light quarks, u, d, and s. A singlet SU(3) symmetry breaking parameter, 
A ~ rri^/m 2 K , incorporates the effects of strangeness suppression in the model, which 
has good predictivity for all the octet 7r and K meson fragmentation functions. The 
model is then extended to include the SU(3) singlet meson, along with additional 
parameters 9p, fd, f g , that describe the singlet-octet mixing. In Section 4, the 
detailed phenomenology is carried out, by fitting the free parameters in comparison 
with known data in both the e + e~ and pp sectors. Section 5 contains a detailed 
discussion and a summary of the conclusions. 

2 Cross sections and rates 

We consider the production of pseudo-scalar mesons in e + e~ annihilation and pp 
scattering processes to leading order in perturbative QCD. We focus primarily on the 
production of (n, K), 77 and 7/ mesons. 

2.1 e + e processes 

To leading order, the hadro-production cross-section in the e + e~ scattering process 
at cm. energy ^/s is [8]: 



1 da h 

(Jtot dx 



J2c q D h q (x,Q 2 ) 
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Here c q are the charge factors associated with a quark qi of flavour i and can be 
expressed [8] in terms of the electromagnetic charge, e^, and the vector and axial 
vector electroweak couplings, Vi = T 3i — 2ejsin 2 6' w and a, = T 3i , as 

47TQ; 2 

= [e 2 q + 2e q v e v q p^s) + (v 2 e + a 2 e )v 2 p 2 (s)} , 

4 = —(v 2 e + a 2 e )a 2 q p 2 (s) , (2) 



Pi(«) 



s(m| - s) 



4 sin # w cos 6 w (m z — s) + m z T z 
1 V 



4 sin cos W / (m 7 — s) + m 7 r 



z 



In Eq. [U T 3 is the third component of weak isospin, 6 W is the Weinberg angle, 
mz, ^ z the mass and width of the Z-boson, and a sum over quarks as well as anti- 
quarks is implied. Here x is the energy fraction, x = E hsiidTon / E heam = 2E h j^fs. 
Since we neglect the hadron masses, this is the same as the momentum variable, 
x p = -Phadron/-Pbcam = 2P h /y/s; x 2 = x 2 + 4ml/ s, for a hadron of mass m h where 
Q = y/s is the energy scale of the interaction. The fragmentation function, D q (x, Q 2 ), 
is the probability at a scale Q for a quark q to hadronise to a hadron h carrying a 
fraction x of the energy of the fragmenting quark. 

We re-express the cross section in terms of the non-singlet and singlet fragmenta- 
tion function combinations, as 

1 da h a Do(x, Q 2 ) + a 3 D 3 (x, Q 2 ) + a 8 Ds(x, Q 2 ) 



a tot dx 



where D , D 3 and D 8 refer to the singlet, D = (D u +Dd+D s ), and the two non-singlet 
combinations, (D u — Dd) and (D u + Dd — 2D S ), respectively, with ao = (c u + q + c s )/3; 
a 3 = (cu — Cd)/2 and as = (c„ + Q — 2c s )/6. Again, a sum over both quark and anti- 
quark flavours is implied. Note that the Dq term dominates at the Z-pole, when 
clq ^> 03, cts- 



2.2 Quark versus gluon jet fragmentation 

The OPAL Collaboration [9J also has data separately on charged hadron produc- 
tion from quark and gluon jets. For the case of charged hadrons in quark jets, the 
normalised rate is given by, 

1 dN c 
N~dx 

where l/o~ tot da h /dx is given in Eq. [3] and the sum over charged hadrons h = 7r + , 
7r~, K + , K~, K s is dominated by the pionic contribution. The factor of 1/2 occurs 
because the rate per quark jet is measured and there are 2 jets per process. 



1 da' 



dx 



(4) 
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The more interesting data is that of charged hadrons from gluon jets. This arises 
in 3-jet production via e + e~ — > qqg where the gluon jet is isolated by tagging on 
the (heavy) quark and anti-quark. The normalised rate here is directly proportional 
to the gluon fragmentation function, since the overall kinematical as well as charged 
factors cancel: 

ldN * =5>f- c« 

9 h 



N dx 



Hence the 3-jet data gives information on the gluon fragmentation function. 



2.3 The pp process 

In addition to the unknown fragmentation functions, hadro-production in pp pro- 
cesses requires information on the initial state parton density distributions. The 
underlying processes here are all possible qq, qg and g g (2 — >• 2) interactions where 
one of the final state partons fragments into the meson of interest. 

The invariant inclusive cross section for the reaction p + p — > h + X for producing 
a hadron h at large pt in the center of mass of the initial state protons (neglecting 
the transverse momentum) is given by [TU] , 

S/iji = - E / dx * / dXb p °( x «> fffa' -^ir- ^ > ( 6 ) 

a Ph 11 Jx™ in Jxf™ Zh.Ot 

where the sum over (a, b, c, d) runs over both quarks and gluons. Here x a and Xb are 
the usual Bjorken-x corresponding to the parent proton momenta pa and ps'- x a = 
Pq/pa, = Pb/PB- Hence P a /A(%a, Q 2 ) are the usual parton density distributions; for 
example, P u / P (x a ,Q 2 ) = u(x a ,Q 2 ), etc. 

The fragmentation functions depend on the variables, z = Zh = Ph/Pc and Q 2 = p\ 
while the partonic sub-process variables are as usual defined in terms of the s, t and 
u hadronic variables t = x a tjz and u = Xbu/z. The limits of integration 

are P2] 

min . min %a%2 

1 X a X a X\ 

with x\ = —u/s,X2 = —t/s. 

For numerical comparison with the data, we reexpress the cross-section in terms 
of the physical observables which are the transverse momentum px = Ph sin 9 and the 
rapidity y = (1/2) \n[(E h + p h cos9)/(E h -p h cos9)], as 



dV 1 d 3 



a 



h dp\ 2pT dpTdyd(p ' ^ 

where 9 is the scattering angle of the hadron h in the pp center of mass frame and 
Eh and ph are its energy and 3-momentum. 

The various sub-process cross-sections [TT] are listed in Table [TJ Note that the 
qq, qg and g g processes all contribute at the same order in a s . Hence the quark and 
gluon fragmentation functions contribute at the same order, unlike in the e + e~ case. 

We now present details of our model for quark fragmentation functions. 
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Subprocess 


eft 2jj 


QiQj -> 'li'lj- 
qiqj -> ft^j 


2F( X ) + 


g<?< -> ft* 


2F( X )-f{ X + 2 + ^} 


-> 

g^gi -> qjqj 
qi^ -> 


2F(x) + {^(x-l + ^) + ^} 

U + TV 1A "I" x (1 +X )2 1 


?s -> w 


/i2F( X ) + {x + 3 + ^} 


^ -> 
-> qq 


^2F( X ) + ^{4-i} 
u+ jv 9 i z U + X J (1+x)2 | 



Table 1: Subprocess cross-sections in pp scattering [TTJ in units of ana 2 /s 2 ; a = 
4Cp/N g , h = C A /C F , with C F = 4/3, C A = 3, iV = 3, and N g = 8. 



3 The Model 

This model was developed in Ref. [7j to study 7r and K fragmentation in e + e~ colli- 
sions. Consider inclusive octet hadro-production under the assumption of exact SU(3) 
(flavour) symmetry: 

qi ^hi + Xj, (8) 

where qi = u,d,s for % = 1, 2, 3. That is, we have 3 — > 8 + X, with X being a triplet, 
antisixplet or fifteenplet. This holds for both octet meson and baryon production. 
The fragmentation of a light quark (u, d, s) into any member of the octet is completely 
described by three independent fragmentation functions, a(x, Q 2 ), /3(x, Q 2 ), j(x, Q 2 ), 
corresponding to X = 3, 6, 15 respectively. 

The corresponding fragmentation functions for the antiquarks are a, (3 and 7 
respectively. As the pseudo-scalar meson-octet contains both the mesons and their 
antiparticles, D^(x,Q 2 ) = D^(x,Q 2 ) so that there are only three independent quark 
fragmentation functions for the entire meson octet, which we choose to be a, (3 and 
7- 

These simplify further on applying equality of unfavoured fragmentation so that 
= Dj + = D~l + = D1 , etc., which reduces the number of independent func- 
tions to a valence fragmentation function, V(x, Q 2 ) and a sea fragmentation function, 
S(x, Q 2 ), where 

V = a + (3--^ = a--^- (9) 
S = 4/5 = 2 7 . (10) 

The detailed expressions for in terms of these functions are given in Ref. p] and 
are reproduced in Table [2] in terms of V and 7. 
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2 7 
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V + 2 7 
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K + 2 7 



Table 2: Quark fragmentation functions into members of the baryon and meson octet 
in terms of the SU(3) functions V(x) and 7(2) in the exact SU(3) symmetric case. 

3.1 Flavour breaking effects 

Flavour SU(3) is badly broken due to the relatively heavier s-quark compared to 
the u- and <i-quarks. This is implemented in the model through an ^-independent 
suppression factor A whenever a strange quark is needed in order to fragment into 
that meson. For example, , + and Df + are suppressed by A compared to the 
SU(3) symmetric expressions given in Table [2] since we require q — > qs in order to 
fragment into K + with a valence s. While and are similarly suppressed, 
only the sea part of is suppressed, so that the fragmentation of s into K + is 
given by (see Table EJ), 

D* + = 2V + 2A 7 . (11) 

The valence component is not suppressed since the heavier strange quark is al- 
ready available and only non-strange quarks are needed to produce K + . The sea 
(anti)-quark fragmentation function is suppressed since, by definition, the fragment- 
ing (anti)-quark is in the sea of the produced hadron and a strange (anti)-quark is 
still needed to make up the valence quantum numbers. Similar arguments can be 
applied to the fragmentation functions of K~, K° and K°. 

In short, the sea contributions to the K meson fragmentation functions remain 
SU(3) symmetric and are uniformly suppressed by A compared to the corresponding 
unbroken SU(3) (or n) fragmentation functions. Hence, all the quark fragmentation 
functions, and hence the combinations Do, D3 and Dg are expressed in terms of the 
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octet functions, V and 7 and the constant suppression factor A, which we expect to 
be A ~ m^/rn 2 K . 

Until now nothing new has been introduced. We now extend this model to include 
the singlet mesons and SU(3) singlet-octet mixing, in order to determine the rj and 
rf fragmentation functions. 

3.2 Extension for 77, rf mesons 

The physical 77 and rf states are orthogonal admixtures of the pure SU(3) octet and 
singlet states, r\ 8 and 771, with mixing angle Op: 

r] = r] s cos 6 P — rji sin Op ; 
rf = rjs sin Op + rji cos Op , 

where \rj 8 ) = (mm + dd — 2ss)/a/6 and \rjx) = (uu + dd + ss)/\/3. The quadratic 
and linear mass formulas predict Op to be —11.5° and —24.6° respectively; present 
experimental limits allow the range —24° < Op < —10° [8]. Note that Op is still small; 
cos Op > 0.9, and so 77 [rf) is dominated by its r] 8 (rji) component. 

The purely octet contributions are given in Table [2J Note that the singlet combi- 
nation of the r] 8 fragmentation function is given by Dq = 2V + I27 just as for the 7r 
and K mesons. The non-singlet fragmentation functions are -D| = and -Df = —2V. 

3.3 SU(3) singlet fragmentation functions 

In order to complete the description for the physical r] meson, we now need to address 
the SU(3) singlet contributions. There is only one singlet fragmentation function that 
we denote as 5, since the only possibility is 3 — > 1 + X: 

qi^h + Xi, (12) 

where qi = u,d,s for i = 1,2,3, with X being a triplet alone. Instead of including an 
additional unknown fragmentation function, we follow a slightly different route. Since 
S describes a pure SU(3) singlet fragmentation function, it should be proportional to 
the singlet fragmentation function D^. We set 

5 = Dl = f d (Dl) = f d (2V+\2 1 ) , (13) 

where we choose fd to be an x-independent (unknown) constant. For the same reason, 
we set the non-singlet functions D\ = D\ = 0. These conditions completely determine 
the individual rji quark fragmentation functions in terms of fd as given in Table [31 

3.4 Mass suppression effects in the 773 fragmentation func- 
tions 

3.4.1 Valence sector 

We now consider SU(3) breaking effects and proceed along the lines discussed for 
suppression in K mesons. Since r/ 8 is a bound state of q i q i quarks of the same flavour, 
only the strange quark valence fragmentation functions, i.e., the valence components 
of -Df and L>| are suppressed; this suppression is again by the same factor A as in the 
kaon case. 
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quark 
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fd&V + Zy) 



Table 3: The unbroken SU(3) singlet rji fragmentation functions; fd is a free param- 
eter. 

3.4.2 Sea sector 

There is an additional complication for the sea quark fragmentation relative to the 
kaon sector: a sea quark has to pick up a quark-anti-quark pair of the same flavour 
so that — > i] 8 , but all possible flavours, i — u,d, s, are possible because of the 
flavour content of the 7]$ meson. If a strange quark-anti-quark pair is produced during 
fragmentation, the contribution is suppressed by A, exactly as in the case of kaons. 

A priori, it appears that there should be no suppression when uu or dd pairs 
fragment into i] 8 . However, a non-strange uu or dd pair will prefer to fragment to 
7r rather than rj because of the relative masses involved. Hence, non-strange sea 
fragmentation to i]g is also suppressed; since the strangeness suppression factor is 
A = m^/m 2 K , we take the non-strange suppression factor to be As = m^/m 2 . Note 
that we have used the relevant mass as that of the physical mass eigenstate, m v , 
since the rj 8 contribution dominates the rj state. Since sea quark contributions to the 
rjs fragmentation can arise from production of quark-ant i- quark pairs of any flavour, 
hence they are uniformly suppressed by a factor / 8 = (2A 8 + A)/3. 

3.5 Mass suppression factors in 771 fragmentation functions 

The logic is exactly the same as for the r/ 8 case. The valence strange sector is 
suppressed by A as before while the sea sector is uniformly suppressed by a fac- 
tor fi = (2Ai + A)/3 where Ai = m^/m 2 , is the non-strange suppression factor for 771, 
assuming that 771 dominates the 7/ state. 

3.6 The 77, 77' fragmentation functions 

Hence, the model has been extended to include a realistic prediction of the r]g and 
771 fragmentation functions in terms of the two octet functions V and 7 and the x- 
independent suppression factor A (all known from fits to 7r and K data), along with 
the ^-independent suppression factors A 8 and Ai apart from an unknown constant fd- 
In summary, we have, for the octet meson: 

D 8 u = D 8 d = ^ + 2/ 87 , 
^ 8 = |vA + 2/ 87 , 
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and for the singlet meson: 

D l u = D\ 

The fragmentation functions for the physical states can then be expressed in terms 
of the fragmentation functions for the rj 8 and rji parts, along with the mixing angle, 
dp. We have 

DV = (c?) 2 (cos 2 O P ^ + sm 2 O P 0^ , (14) 

D i = (4? (sm 2 Op^I + cos 2 , (15) 

where i = u,d,s and the coefficients are = c d = 1, c® = —2, c l u = c l d = c\ = \^2, 
c u = c d = (cos6p — \/2sm8p), and c v s = (— 2 cos#p — v^sin^p). The coefficients for rj 
are obtained from those of r\ by the substitution: (cos 9p — > sin Op] sin 9p — > — cos 9p). 

4 Comparison with the data 

We now fit the unknown model parameters by a comparison with data. Note that the 
octet fragmentation functions V(x, Q 2 ) and 7(x,Q 2 ), the suppression factor, A, and 
the gluon fragmentation function D g (x, Q 2 ), have already been fitted to the 7r and K 
data in Ref. [7]. We consider them here again for reasons explained below. 

4.1 e + e — > 2-jets 

The V and 7 were fitted to parametrised forms as a function of x: 

Fi(x) = Oi(l - x) b *(x Ci )(l + diX + eiX 2 ) , (16) 

at a starting scale of Qq = 2 GeV 2 and evolved, along with the gluon fragmentation 
function, D g (x, Q 2 ), to leading order (LO) in perturbative QCD. Best fits to V(x, Qq), 
S(x,Qq), D g (x,Ql) and A were obtained from fits to electro-production data at the 
Z-pole Q 2 = 91.2 2 GeV 2 for tt ± and K± [12]. 

The best-fit values for the parameters a, b, c, d, e for different input fragmentation 
functions are given in Table HI The best fit value for A is 

A = 0.08 , (17) 

which is consistent with the quark model expectation A = m^/rri^. 

The resulting fits to the ix and K production rates at the Z-pole from LEP data 
[I~3| IT4] are shown in Fig. [lj Note that the gluon is not well-determined from this 
data, unlike the quark fragmentation functions. This is because the cross-sections are 
functions only of the quark fragmentation functions at LO and the gluon contribution 
only enters through the evolution equations. The electro-production rates are insen- 
sitive to the gluon even at next-to-leading order (NLO) since the gluon contribution 



/ d ( -y + 2/ l7 



f d [-VX + 2f 11 
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is down by a factor of a s compared to the quark contributions. Three-jet processes 
are sensitive to the gluon contribution; we return to this point later. 

Note that the sea quark and gluon distributions are different from those obtained 
in Ref. [7j. Since gluons and sea quarks mix during singlet evolution, a small non- 
zero gluon was included earlier while evolving from the starting scale Ql = 2 GeV 2 to 
the observed scale Q 2 . Gluon fragmentation function data is now available [9] from 
isolating charged hadrons produced from gluon jet fragmentation in e + er — > qqg 
processes. It turns out that the gluon fragmentation used in Ref. [7J is too small and 
therefore incompatible with the 3-jet data. Our new updated fits here include a gluon 
that is consistent with such data (as discussed below); this necessitates a revision of 
the sea quark starting fragmentation functions (they have been made smaller) so 
that the sea contribution evolved to the Z-pole (with the contribution from the new 
enhanced gluon) still matches the earlier result (and hence the data). 

The valence function V{x, Qq) remains the same as before. 




lO- 2 10- 1 lO- 2 10- 1 1 

X X 



Figure 1: The model fits to the Z-pole data for 7r° [13] and K° + K [T3] mesons; the 
functions at the starting scale Qq = 2 GeV 2 are given in Table HI 

It is observed that the predictions exceed the data systematically at small-x. This 
is because of the pole in the P gq and P gg splitting functions as x — > which cause 
the DGLAP evolution to drive the small- x sea and gluon fragmentation functions 
to large valueswhich cause the DGLAP evolution to drive the small- x sea and gluon 
fragmentation functions to large values. While the starting values of the fragmenta- 
tion functions are integrable, this evolution makes them non-integrable due to this 
increase at small-x. However, the fragmentation functions should remain integrable 
(as their first moments are related to the hadron multiplicities). An approach such 
as the modified leading log approximation (MLLA) [15] cures the small-x divergence 
and gives a good fit to the data at small- x as has been discussed in Ref. [7]; here 
we simply ignore the small-x region and concentrate on the fits for x > 0.02 for all 
mesons under consideration. 
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2.33 


3.5 


7.25 


b 


2.15 


12.76 


4.4 
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-0.64 


-0.75 


-0.5 
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5.35 


3.87 
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e 


-5.12 


61.59 






Table 4: Input values at Q 2 = 2 GeV 2 for the valence, sea and gluon fragmentation 
functions for the pseudo-scalar meson octet. The factor of 2 is due to the convention 
used in the earlier work [7] when the model was fitted to the sum of the charged 
hadrons, for example, ir + + 7r~, K + + K~ , etc. 

4.1.1 Fits to r], rj' data on the Z pole 

With As and Ai fixed to the theoretically expected values, fd is the only unknown 
parameter in the rates for both rj and rf production. Since the nonet mixing angle 
6p is small, fd is mostly constrained by the large-x data. Fits to the rj, rj' production 
rates in e + -e~ annihilation on the Z-pole at LEP [T6], [17] yield a best fit value to the 
free parameter of fd = 0.3. However, the electro-production data are insensitive to 
the nonet mixing angle Op, as can be seen from the fits to the data in Fig. [2] where the 
fits for 9p = 0°, —15° and —24° are shown. The allowed range of 9p marginally alters 
the fits at small-x values, while being completely insensitive in the large-x region. 




Figure 2: The figure shows the model fits at y/s = 91.2 GeV to the LEP data on rj 
[j~6] and [17] mesons as a function of x with both x and y error-bars. The green, 
red and blue lines represent Op = 0°, —15° and —24° respectively for fd = 0.3. 

4.2 pp — »■ (it, rj)X 

The PHENIX experiment at RHIC [21 [1] has measured the inclusive transverse mo- 
mentum spectra of tt and r] mesons in the range pt ~ 1.5-14 GeV/c at mid-rapidity 
(\r}\ < 0.35) in pp, dAu and AuAu collisions at ^/snn = 200 GeV. They have also 
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measured the 77/71"° production ratio, R^/^o. This ratio is interesting because of its 
potential as a good signal for quark-gluon-plasma. 

The 77/71"° ratio in pp collisions is the base-line with which comparisons in d A 
or A A processes are made. It is therefore important to understand and correctly 
determine this ratio m pp collisions for a wide range of Q 2 (~ p 2 ^). As can be seen 
from Eq. [61 this depends on the initial parton density distributions and the final state 
fragmentation functions. Furthermore, the gluon and quark fragmentation functions 
(and density distribution functions) occur at the same order in a s . The study of 
such a process, therefore, may help determine the gluon fragmentation functions with 
much greater accuracy than currently known. With these two considerations in view, 
the model was compared to the rj and 7r production rates in pp collisions. 



4.2.1 Analytical approximation 

Before comparing the model to the data, it is observed that there is considerable sim- 
plification of the expressions when the pp sub-process cross-sections are reexpressed 
as a constant times a kinematic factor, F(x), with small correction terms [TT], as can 
be seen from Table [TJ, where 

F( X ) = X 2 + X + l + - + ^, (18) 

X X 

where \ = u/t. When these correction terms are dropped, the surviving terms are 
such that the fragmenting parton c has the same flavour as the initial parton a. Then, 
apart from an overall factor, 

the integrand of Eq. [6] for the hadro-production rate can be factorised as the product 
of two terms: 

(<hD} + QiD- + h gD^j x ( qj + q 3 + hg) . (20) 

i,3 

The first term depends on (x a , z) and contains the dependence on the unknown frag- 
mentation functions while the second term depends on Xb alone. In this factorised 
form it is clear that there is a substantial dependence on the gluon fragmentation 
function only when the corresponding gluon density distribution function is large as 
well; this occurs, for the central rapidity region measured by PHENIX at RHIC, for 
low p T , p T < 5 GeV. Hence, the small-p-r data (with Q 2 « p\ > 2 GeV 2 ) should be 
sensitive to the gluon fragmentation function. Since our input fragmentation func- 
tions are at a scale Qq = 2 GeV 2 , there may be significant threshold effects at this 
energy; we choose Q 2 = (p? + ml) with Q 2 > Qq. 



4.3 Fits to the tt production rate 

The expressions for the production rate are integrated over the central rapidity region, 
—0.35 < y < 0.35, and over the allowed ranges of x a and x^ for a given p?. The fits 
to the 7r production data [2] as a function of pt are shown in Fig. [3J Here the exact 
expressions for the sub-process cross-sections have been used, not the approximations 
discussed in the previous sub-section. The small-pr data, as expected (where Zh > 
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0.02), are very sensitive to the gluon fragmentation function, while the larger p p data 
( z h ^ 0-1) are sensitive to the quark fragmentation functions. The effect of scale, 
Q 2 = fsipr + m 1)i fs — 0.5-2, is also shown in the figure. 




123456789 101112 
Pt 

Figure 3: The 7r production rate in pp collisions at y'i = 200 GeV in comparison 
with the RHIC/PHENIX data |2J. The central curve corresponds to Q 2 = p\ + Tn\ 
while the upper (lower) curve corresponds to a value of Q 2 half (twice) this value. 



4.4 Fits to the r] production rate 

It is reasonable to assume that the gluon fragmentation function in rj is suppressed 
relative to that of tt, just as is the case with sea quarks. However, there is a difference 
during evolution: the gluon mixes with the entire singlet combination, Dq, which 
contains both valence and sea quark terms. For example, 



where 



Dl = 2K+12 7 
D k = (i±^) 2 y + A12 7 
Dq = x 2\/ + y 127 



x = I [2(cl) 2 + A(c^) 2 ] (cos 2 6 P + sin 2 6 P ) , 

ep ^ lX -0.5; (21) 

Vo = ^ (cos 2 6 P [8(cD 2 + (c^) 2 ] + f d f lS m 2 9 P [4(c^) 2 + 2(c^) 2 ] ) , 
e P small , 

- * J8 • 
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Hence the valence in D ,r? is suppressed by a factor of 1/2 and the sea is suppressed 
by A, f 8 compared to Dq . Therefore, we expect the gluon fragmentation into r? to be 
suppressed by D v g ~ Df = f g D g , where y < f g < x . 

Since V and 7 dominate in different regions of phase-space (at different Zh values), 
it is likely that f g is an x-dependent suppression factor. This can be studied in 3-jet 
processes in e + e~ collisions (see the next sub-section) . Since detailed data on this is 
still not available, we assume f g to be an x-independent constant between y and x ; 
certainly f g < 1. 

The resulting fits to the rj production rate in pp collisions are shown in Fig. H]for 
different allowed mixing angles Op. While the data is sensitive to f g , as expected, 
it does not appear to be very sensitive to Op. Note that the maximum separation 
between the curves for a given mixing angle is at small-p^, where the gluon con- 
tributes maximally. At large enough pp, the ratio should determine Op as the gluon 
contribution becomes negligible. 




Figure 4: The 77 production rate in pp collisions at y/s = 200 GeV as a function of 
p T with Q 2 = p\ + m\ in comparison with the RHIC/PHENIX data pQ. 

4.5 Fits to the 77/71"° ratio 

The sensitivity to both f g and 9p is amplified in the t]/tc ratio, as can be seen in 
Fig. [5] where the RHIC/PHENIX data [1] are also superposed. The region bounded 
by the curves for f g = y ,lin the figure encompass the possible fits to the ratio. 

A x 2 minimisation gives a best fit for f g = 0.3, Op = —20.0°. Since the total errors 
are large, only the statistical errors were used in the fit to give x 2 = 9-5 for 11 d.o.f. 
A fit with the total errors reduces this best fit value to x 2 = 4.5. The saturation 
value of the r]/7r ratio increases with increasing \9p\. Constraining the large pp ratio 
to be ~ 0.5 p] as favoured by PYTHIA [18] decreases the value of Op to about —16°; 
however, the errors on the data also allow this ratio to saturate at values of around 
0.6. An improvement in the data quality will definitely constrain both Op and f g 
better. Currently, a reasonable quality of fit (x 2 ~ d.o.f.) is obtained in a range 
f g = 0.3-0.35 and —24° < Op < —16°. Also, while the 7r and rj production rates are 
scale dependent, this dependence is very small and virtually cancels out in the ratio, 
except for a mild dependence at smaller pp values, pp < 3 GeV. 
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Figure 5: The 7] /it production ratios in pp collisions as a function of pt for different 
nonet mixing angle Op for two different gluon suppression factors, f g = yo,l, in 
comparison with the RHIC/PHENIX data at y/s = 200 GeV pQ 



The quality of fits to the rj/n ratio for the typical choices: (f g , 6 P ) = (0.3, —20°) 
and (0.35, -16.77°), are shown in Fig. © along with the RHIC/PHENIX data where 
both statistical and total errors are shown. 
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Figure 6: The f]/iT production ratios in pp collisions as a function of px for 
(f g ,0) = (0.3, —20°) (solid line) and (0.35, —16.77°) (dashed line) in comparison with 
the RHIC/PHENIX data [lj at ^ = 200 GeV. 



4.6 Fits to the 3-jet data in e + e collisions 

Fits to the individual quark fragmentation functions are of the same quality as for 
the 7r, K production discussed earlier (this data is just the sum of the 7r and K rates). 
Hence we do not discuss them further and proceed directly to gluon jet fragmentation. 
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With f g = 0.3 as determined from fits to the rj/ir ratio in pp collisions, it is now 
possible to compute the charged hadron production rate from gluon fragmentation in 
e + e~ — ► qqg. We have, from Eq. [5j 



1 diV ch 



N dx 



2.9-3.2D* 



where we have assumed that fy- = fJJ = 0.3-0.35, so that the gluon fragmentation 
into r) and K are roughly suppressed by the same factor. If ir and K production from 
gluon jets is separately measured, f g (and its x- dependence, if any) can be easily 
determined. Such data would be very useful in constraining further the small-pr 
behaviour of the r]/n ratio. 

The resulting fit to the OPAL data [9] is shown in Fig. [7J Notice that the Q 2 scales 
of the e + e~ and pp processes are very different: the former with a value Q 2 ~ 80 2 
GeV 2 (slightly smaller than the Z-pole value for 2-jet production), and the latter 
with Q 2 = p\ = 2 2 -14 2 GeV 2 . Due to the pole in the P gg splitting function, the 
gluon fragmentation function is very different (in fact, at both large and small Zh 
values) in the two cases and gives a strong consistency check on the fits to the gluon 
fragmentation function. 



X 




Figure 7: The normalised charged hadron production rate from gluon jets in e + e~ — > 
qqg at a mean Q 2 = 80 2 GeV 2 at LEP/OPAL [9] along with the theoretical prediction 
Dg ~ 2.9-3.2Dg where h runs over all charged hadrons and Ks- The lower (upper) 
curve corresponds to the choice of 2.9 (3.2), although the separation between them is 
not well-marked. 
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5 Discussion and Conclusions 



Quark and gluon fragmentation functions of n, K, rj, and rj' mesons have been de- 
termined through a study of hadro-production rates in both e + e~ and pp collisions. 
Detailed information on the r\ fragmentation functions has become important since 
the variation of the rj/n production ratios in pp, pA and A A collisions is a possible 
signal of quark gluon plasma (QGP). 

All calculations have been performed at leading order (LO) in QCD. Inclusion of 
NLO terms should not affect the quality of the fits, especially that of the gluons, since 
they have been mostly determined from ratios (normalised charged hadro-production 
in e + e~ — > 3-jets and tj/tt production ratio in pp collisions). Of course, numerical 
details, such as the fits to the quark and gluon fragmentation functions at the starting 
scale Qq = 2 GeV 2 , will change. The fits therefore are mainly driven, and limited by, 
the model assumptions. 

A simple SU(3) model for octet meson production (jr,K), has been extended 
to include 77 and rj' production. Due to SU(3) symmetry, the quark fragmentation 
functions of all mesons are related to a common valence V(x, Q 2 ) and sea quark 
7(2;, Q 2 ) fragmentation function. SU(3) breaking due to the heavier strange quark is 
included through a constant parameter A. Inclusion of singlet-octet mixing causes 
SU(3) breaking in the sea sector as well; this is parametrised by a nonet mixing angle 
9p that is known to be small (—24° < Op < —10°), and a few additional constants, 
A 8 , Ai, and f d . 

While A, A 8 and Ai are fixed from model considerations, the data, especially on 
the Kj 7r ratio in e + e~ collisions, also yields the same value of A as that predicted 
by the model. The parameter fd that occurs in the singlet fragmentation functions is 
fitted from rj and rj' data in e + e~ collisions on the Z-pole. 

Finally, suppression of gluon fragmentation in K, rj relative to 7r is parametrised 
by f g (x). In the absence of more detailed data (for instance, on individual n, K 
production from gluon fragmentation in 3-jet processes at LEP), f g is taken to be 
constant and fitted from the large-pr (2 < p? < 11 GeV) rj/n ratio in pp collisions 
at y/7s = 200 GeV. 

This gluon parametrisation, with f g ~ 0.3-0.35, is completely consistent with 3- 
jet data in e + e~ collisions at a very different scale Q 2 ~ 80 2 GeV 2 , where the gluon 
fragmentation is explicitly studied. 

A comprehensive set of gluon and quark fragmentation functions (for all favlours) 
is therefore available for the set of nonet pseudo-scalar mesons, n^' , K^, K°, K , rj, 
and rj', starting with a simple model with reasonable assumptions, that is consistent 
with hadro-production data in both e + e~ and pp collisions. The focus of this paper 
was the extension of the earlier fits [7] in 7r and K to the case of 77 and rj' mesons. 

Acknowledgement : We thank H.S. Mani and M.V.N. Murthy for many detailed 
discussions. 
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